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1. Introduction 

Theoretical investigations of spin systems are often motivated by the apphcations to real 
systems as magnetic molecules or spin lattices. In the case of a system of three spins 
with realistic spin quantum numbers varying from s = | to, say, s = |, the underlying 
Hilbert space has a maximal dimension of (2s + 1)'^ = 8'^ = 512. Hence, even without 
using any dimensional reduction due to rotational symmetry, it is an easy exercise to 
numerically calculate physically relevant quantities, e. g. , specific heat or magnetization 
curves, by using standard codes. Other aspects of spin triangles have been treated in 
recent publications [1], [2], and [3]. Thus a detailed theoretical investigation of the gen- 
eral spin triangle requires additional justification. In my opinion, the reasons to study 
this special class of spin systems are threefold. 

First, the spin triangle is interesting in its own right since it exhibits salient fea- 
tures which might evade the purely numerical treatment. It is a straight forward task 
to reduce the pertinent Hamiltonian H to that of a two-chain with coupling constants 1 
and y and to restrict it to invariant subspaces of maximal dimension 2s + 1, see section 

2. Here the reduced Hamiltonian has a tridiagonal matrix representation, see section 3, 
which can be understood by virtue of the Wigner-Eckhard theorem. Some properties of 
tridiagonal matrices that are relevant for our purpose are collected in Appendix A. The 
corresponding characteristic polynomial and the spectral curves e^{y) exhibit a number 
of symmetries that I have dubbed "spectral symmetries" in section 4. Usually, symme- 
tries in quantum mechanics lead to operators commuting with H. Spectral symmetries 
are of a different kind and can be expressed by functional equations for the ey{y). They 
can be traced back to permutational invariance, not of a particular Hamiltonian, but of 
the whole class of triangle Hamiltonians. Despite its elementary origin, spectral sym- 
metry allows a qualitative discussion of the spectral curves and, e. g. , the calculation 
of the moments of the restricted Hamiltonian up to fifth order without using its explicit 
matrix representation. 

A second reason to study spin triangles is that their properties might help to under- 
stand more complex systems. One example, not considered in this article, is the = 
ground state (for integer s and approximately equal coupling constants) that has been 
used to build ground states of larger spin systems including spin lattices, see [4] and [5]. 
Another example is the observation that the first few terms of the perturbation series 
for the eigenvalues of H yield certain classes of universal polynomials that occur in the 
graph expansion of the /c-th moments of if, see section 5 and 6. These polynomials 
also occur in the high temperature expansion for general Heisenberg systems, not only 
triangles, see [8] and [9] . Appendix B contains explicit expressions for the fc-th moments 
of spin dimers that are used in the formulas for the the above-mentioned polynomials. 

Finally, spin triangles may serve as simple but non-trivial examples to study 
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theoretical problems of general interest. Here we will concentrate on one aspect, namely 
integrability. Classical systems are completely integrable in the sense of Arnold-Liouville 
iff they possess Poisson-commuting integrals of motion, where 2N is the phase space 
dimension. Sometimes the corresponding quantum systems are also analytically solvable 
(harmonic oscillator, hydrogen atom) but not always. One counter-example is the one- 
dimensional particle in a potential V{x), another one is the general spin triangle. In the 
latter case we have classically three unit spin vectors, i. e. , a six- dimensional phase 
space, and three constants of motion, H, S"^ and 83. The corresponding quantum 
system, notwithstanding the partial results of this paper, cannot be solved in closed 
form. However, there are weaker notions of "quantum integrability" that justify, e. g. , 
to consider the s = ^ chain as integrable by means of the Bethe ansatz , see [6] or 
[7] for generalizations. In this context one could see the results of the present paper 
as another example within the wide range of the not yet fully understood notion of 
quantum integrability. We will come back to this question in the section 7 containing 
the summary and an outlook. 



2. Basic definitions and elementary results 

We consider the general Heisenberg spin triangle with one and the same individual spin 
quantum number s = |,l,|,2,... and the Hamiltonian 

H = J1S2 ■ S3 + J2S3 • si + J3S1 ■ S2 (1) 

where the Ji, J2, J3 are three arbitrary real parameters and the Sj denote the three 
individual spin vector operators, i = 1,2,3. The Hamiltonian (1) can be analytically 
diagonalized if two of the Jj's coincide (i. e. for the isosceles triangle) or if s is small in 
a sense to be made more precise subsequently. In the first case of, say, Ji = J2 = J we 
have 

H = Hj = J{si + S2) ■ S3 + J3 si ■ S2 (2) 



where 

S = si + S2 + S3 , (4) 

5i2 = si + S2 . (5) 

It is well-known that S"^ and commute and have simultaneous eigenvalues of the 
form S{S + 1), S'i2(5'i2 + 1) where 

S = 3s, 3s — 1, . . . ^ (6) 

5i2 = 2s,2s-l,...,0, (7) 

and the triangle inequality 

\Si2-s\<S<S,2 + s (8) 
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holds. From this the eigenvalues of Hj immediately follow from (3). The corresponding 
eigenbasis of Hj will be used to split the matrix representation of the general H into 
smaller block matrices. 
In the general case we write 



^= Y(5^-3s(s+l)j+(J3-J2)(^Si-^2 + -^— ^S2-S3j , (9) 

and, since H and S"^ commute, the problem of diagonalizing H can be reduced to the 
task of finding the eigenvalues and eigenvectors of 

H2{y) = si- S2 + ys2- S3, where (10) 

Jl — J2 . N 

y =-j — J- (11) 

in the invariant subspace 'H(S', 6*3). The latter denotes the common eigenspace of 
and 5*3 corresponding to the eigenvalues S{S + 1) and, say, 

^3 = ^ • = I ^ ■ ^ i^^^S^^ (12) 
™° 1 i • half-integer 

It will be simply denoted by 7/(5'), since the eigenvalues of -^2(2/) do not depend on the 
chosen eigenspace of S3. The dimension of 'H(S') can be determined by counting the 
possible values of the quantum number S12. We obtain 

s - S < S12 < S + s for 0<S<s, (13) 

S-s<Si2<2s for s<S< 3s, (14) 

according to the triangle inequality (8). In order to avoid the case distinction (13), (14) 
in what follows, we will write the above condition in a uniform way as 

9i < 9 < 92 , (15) 

where g = S12 and 9i{s, S), g2{s, S) are defined according to (13) and (14). 
The dimension of 'HiS) increases in steps of two in the range < 5 < s and decreases 
in steps of one in the range s < S < 3s. The maximal dimension A'max = 2s + 1 is 
assumed for S = s. The figures 1 and 2 show typical examples. Especially, it follows that 
for s < 3/2 the dimension of 'H(S') will be at most four and hence H can be analytically 
diagonalized. At this point one should mention that for moderate s, also exceeding 
s = 3/2, the computer-algebra system MATHEMATICA renders the eigenvalues of 
H2{y) in the form of root objects which have many properties of special functions. 
For example, the n-th eigenvalue of H2{y) can be expanded into an exact power series 
w. r. t. y up to some finite order. 

The total dimension dtot of the eigenspace of 5*3 with eigenvalue S'l-^jn amounts to 



c^tot = 3s(s + 1) + <( 3 _ , (16) 



1 : s integer 
^ . s half-integer 

in accordance with the general formula in [10] 
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s=2 



Figure 1. The range of the quantum numbers 6*12 and 5 for s = 2. It is typical for integer s that 
S has a minimal value of S'jj^jj^ = and the corresponding eigenspace 'H{S) will be one-dimensional. 



Note that dimensions of the various spaces 'H(S') are 1 + 3 + 5 
accordance with (16). 



4 + 3 + 2 + 1 = 19 = 3s(.s + 1) + 1 in 
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Figure 2. The range of the quantum numbers S12 and S" for s = 5/2. It is typical for half-integer s that 
S has a minimal value of S'jjjjj^ = 1/2 and the corresponding eigenspace 'H{S) will be two-dimensional. 



Note that dimensions of the various spaces 'H{S) arc 2 + 4 - 
in accordance with (16). 



-6 + 5 + 4 + 3 + 2 + 1= 27 = 3s(s + l) + 3/4 



3. Matrix representation of the reduced Hamiltonian 

In order to find the matrix representation of H2{y) within the invariant subspace 
'H(S') one has to choose an appropriate basis. An obvious choice is the eigenbasis of 
^'2(0) = ■ S2 = \ (yS\2 — 2s(s -|- The transformation to the eigenbasis of S2 ■ ^3 
within HiS) is achieved by the unitary, real, symmetric matrix U with the entries 



(17) 



where (. . .) denotes the 6 J— symbol, see, e. g. , [11]. For the sake of simplicity we 
will not distinguish between operators and their matrix representation in 'H(S'), if no 
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misunderstanding is likely to arise. Hence we may write 



si-S2 = diag{-(^((7 + 1) - 2s(s + 1))|^ = gi, ■ ■ ■ , g2} (18) 
and 

S2 ■ S3 = f/* Si ■ S2 f/ . (19) 

Usually, the matrix U is fully occupied and contains only few zeroes. Thus it is somewhat 
surprising that S2'*3 will always be tridiagonal. On the other hand it is a straight forward 
task to directly calculate S2 ■ S3 in the eigenbasis of Si • S2 if one adopts the theory of 
irreducible tensor operators (ITO), see, e. g. [12] or, for a recent account, [13]. In this 
theory the tridiagonal form of S2 ■ S3 is a direct consequence of the Wigner-Eckhard 
theorem, since S2 is an ITO of rank one. Moreover, it follows by the same theorem that, 
for S = s, (s2 ■ S3)oo = which we will derive later by different arguments. We will not 
dwell upon the details of the ITO calculations and simply present the result. In order to 
avoid the case distinction (13), (14) in the formula for the secondary diagonals of H2{y) 
we will express s and 5" through gi and (72, the bounds of g, see (15): 

^'11 + 92), S = l{g2-gi), for 0<S<s, (20) 




for 



for 0<S<s, 
s<S<3s. 

y 



2s{s + l)) + ^{-g{l + g) 



s{l + s) + S{S+l)) 



9i){9 - 92){g - gi - g2)i2 + g + g2){2 + g + gi + g2) 
16{l + 2g){3 + 2g) 

The other matrix elements of H2{y) vanish. For example, if s 
the matrix form 
/ -12 A^/3y 



4^31/ 










-11 - 

2 







3V15^< 
2 

^ 2 












Q^fly 

-6 - ?>y 

lly 










lly 



5y 



10y_ 

Vn 









VTi 

I5y 



(21) 

(22) 
(23) 

(24) 

S = 3 then -^2(2/) has 
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(25) 



We note in passing that a negative sign of the square root in (24) gives the matrix rep- 
resentation of H2{y) = si-S2 + y Si-S3. This will be later used in section 6 to derive (90). 



The diagonalization of tridiagonal matrices is considerably simpler than for general 
matrices and has been extensively studied also from a numerical point of view, see. 
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e. g. [14]. For the related problem of the inversion of tridiagonal matrices see [15]. We 
have collected some results which will be of relevance for our purpose in the Appendix 
A. For example, in our problem the entries (24) of the secondary diagonals do not 
vanish. In this case the eigenvalues of the tridiagonal matrix will never be degenerate. 
Of course, there is the usual degeneracy due to rotational invariance and some accidental 
degeneracy between the energy eigenvalues of H of different invariant subspaces "^(5*), 
not to mention the special case of Ji = J2 = Js- 

Another interesting property is the following : If we perturb a diagonal matrix by a 
tridiagonal perturbation the resulting perturbation series will be given by relatively 
simple formulas, at least for the first few terms, see (114)-(116). These formulas seem to 
be of little use at first sight since, at least for small s, the eigenvalues can be numerically 
calculated very fast, but we will present different applications in the sections 5 and 6. 



4. Spectral symmetries 

For given quantum numbers s and S let p{x, y) denote the characteristic polynomial of 
H2{y) restricted to 'H(S'), that is, 

p{x, y) = det (^2(y) -X 1 ms)) . (26) 

It is a polynomial in the variables x and y of total degree = dim('H(S')). For example, 
ifs = S' = 3/2we have A^ = 4 and 
135 

p{x, y)= - —{y - 1)^(33 + 94y + 33^/2) 

-:^(|/ + l)(47-110y + 47y2)x 
lo 

_^(l_26y + yV 

o 

+ 5(1 + y)x^ 

+ (27) 

Obviously, the above coefficients of y'^ x^ are symmetric under the reflection v ^ 
N — fi — u or 

N I X 1 



r>(-,-) = p{x,y)ioTy^O. (28) 
This can be proven generally as follows. 

The two Hamiltonians -^2(1/) = ' £2 +2/ S2 ■ S3 and y H2 (jj^ ) = y ^si -82 + ^ S2 ■ S3 j 
are unitarily equivalent under the permutation (13) and hence have the same eigenvalues, 
say, eyiy), u = 1, . . . , N. This means that, upon an appropriate ordering of eigenvalues, 

e.(2/) = 2/e, , (29) 
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from which (28) follows by 



N 



= nf^e-f^V^) (31) 



N 

.N 



u=l 



y n(^^ ( 7 ) - 7) (32) 



yj y 



= y%{-.-) . (33) 
\y y) 

The coefficients of %p are determined by the well-known eigenvalues of -^2(0) 

e,(0) = ^(^7(^7 + l)-2s(s + l)), g = g^,...,g2. (34) 
For example, the factors in front of the lines of (27) result from 

XI I — ■ x\ |— - — XI I - — X 



4 yv4 /v4 yv4 

135-33 9-47 9 2 r 3 4 

= X X +5x -|-x . (35) 

256 16 8 ^ ' 

There exist further symmetries of p(x, y) which are not so obvious as the reflection 

symmetry. They are revealed by employing arbitrary permutations of the spin triangle. 

These permutations can be composed of transpositions [if) that induce rational 

transformations oiy = j^ZjI the following way: 

(13) 1 (12) 1 y-1 (23) y (13) y 1 (12) (23) . 

y — y ^ 1 - - = > — y 1 = yl-y — yy ■ (36 

y y y y-i 

For almost all real numbers y the permutation group S3 thus produces cycles of six 
different numbers. The exceptions are 1 — t- 1 — — )■ ±00 — )■ ±00 — )■ — > 1 and 
|— ;'2— )■— 1— ;■— 1— )-2— )-|— )-|, where only three different numbers are involved. In 
general, if y runs through such a 6-cycle the spectrum of the corresponding Hamiltonians 
H2{y) will be related by affine transformations. (29) is a first example. A second 
example which, together with (29) generates all remaining ones, is 

e,(l -y)=Eo- e,{y) , (37) 

where 

Eo = ^(5(5+l)-3s(s + l)) , (38) 
and a suitable ordering of the eigenvalues is assumed. 

Geometrically, (37) means that the system of spectral curves y h- )■ ey{y), v = 1, . . . ,N 
in the y — -E-plane is point-symmetric w. r. t. the point P = [y = ^, E = ^Eq). For odd 
N, the curve y 1— )■ eN+i{y) passes through P and is point-symmetric in itself. Especially, 
E = ^Eq is always an eigenvalue of -^2(|) for odd and the other eigenvalues of -^2(|) 
lie symmetric to the value ^Eq for all A^. These properties are illustrated in figure 3. 
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E 




Figure 3. Spectral curves ey(j/), = 0, . . . , 8 of a 2-chain with coupling coefficients Ji = 1 and J2 = V 
and s = iS = 4. The system of spectral curves is point-symmetric w. r. t. the red point P = (1/2, —10). 
For pairs of y- values connected by the transformation (36) the corresponding eigenvalues ey{y) are 
related by an affine transformation. This is illustrated here for the pair yi = 1/2, 1/2 = 1 — = ^1 
by demonstrating that the (dashed) lines joining the corresponding eigenvalues eu{y) meet at the point 
(2,-20). 




Figure 4. The same system of spectral curves as in figure 3, but with a larger y-domain. For large 
\y\ the spectral curves approach certain (dashed) lines which are determined by the first two terms 
of the perturbation series of 6^(1/) via (43). For y — +00 the slope of the asymptotic line equals 
the unperturbed eigenvalue (0); for y — >• —00 the slope is e7v_2y(0). For y > the convergence of 
the spectral curves to their asymptotic lines becomes slower for the lower curves due to larger values 
of lajj"^^]. The lowest spectral curve with a large (negative) curvature at y = shows the slowest 
convergence. For y < analogous remarks applies to the upper curves. 
The asymptotic lines meet at the red point P = (1/2, —10). 
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A related consequence of (37) is that the shifted characteristic polynomial 

q{x,y)=p(^x+Y^y + \^=Yl<l'^^y^^'' (39) 

will contain only even or odd sub-polynomials Quiv)- For example, the shifted polynomial 
pertaining to (27) reads: 

g(x, y) = 7^(945 - 5130 - 4455 y') -^{27y + 423 y^) 
25o lo 

- ^{63 + 9y'^) + 5x^y + x^ . (40) 
8 

Unfortunately, these and the above-mentioned symmetries (29) are not sufficient to de- 
termine the coefficients of p{x, y) completely. 

We will explore some further consequences of (29) and (37). For \y\ — j- oo the 
spectral curves y H- Cuiii) will approach certain lines. This can be easily understood by 
inserting the perturbation series 

e,{z) = x'^^^ + x^^^ z + x^^^ + . . . (41) 
into the r. h. s. of (29) and setting z = ^ 

V y y J 

= x^Q^ y + x^^^ + x^^^ ^ + . . . , for y ^ +00 . (43) 
For y — )■ — oo one obtains 

e.(y) = x^-^^^ y + x!^-^) + xf - + .... (44) 

y 

Since the spectral curves do not intersect (see Appendix A), it follows that for y — )■ +oo 
the spectral curve ey{y) approaches the line with slope x^q^ and for y +oo it 
approaches the line with slope x'"^~'^\ This effect is illustrated in figure 4. 
Anticipating the result (50) it can be shown that the asymptotic lines meet at the 
point P=(2/ = i,E = iEo): 



By the sequence of rational maps (36) the closed interval [0, |] is mapped as follows: 
[0, ^] ^ [2, +oo] ^ [-00, -1] ^ [-1, 0] ^ [1, 2] ^ [i, 1] ^ [0, ^] . (46) 

This set of six closed intervals covers the whole real axis, hence the values of e^{y) 
for y G [0, |] completely determine the full spectral curve e^{y), ?/ G M. We have 
already exploited this principle for the discussion of the asymptotic behavior of e^{y) 
for \y\ — oo. 

Conversely, if we choose some analytic function, restrict it to ?/ G [0, |] and extend it 
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to the full real line by means of the functional equations (29) and (37), the resulting 
function will not automatically be analytic at the borders of the intervals, e. g. at y = 0. 
This makes it plausible that (29) and (37) impose conditions on the perturbation series 
(41). Let us consider some small y. After the first three steps of the sequence (36) y is 
mapped onto which is again small. Simultaneously, e^i^y) is mapped onto 

Eo-, (i^^o - e^{y)) = - E^-^ + Y.'^l'V (47) 

1 — V 1 — V 1 — V ^-^ 



=0 



Equating the coefficients of the power series (47) and (48) leads to the following system 
of equations: 

= -Eo + f2(^ + i-^r'('^Zl))^n\ m=l,2,... (49) 

n=0 ^ \^ / / 

Consequently, for all u = 1,...,7V the perturbation coefficients Xn \ n odd, can be 
expressed in terms of the perturbation coefficients xf\ i even and i < n. The first few 
of the corresponding equations read: 

I {Eo - ) (50) 

(51) 



X 



2^5 



X 

W 

X. 



-x''^^ + 6x1"^) (52) 
xt'^ - 5x1") + ^^t^) (53) 



1 

4 
1 

2 

17 (^) 21 (^) 35 (^) 7 (^) 
- y4 + - -^A + • (54) 

Sometimes we will switch from the index notation of eigenvalues using v = 1, ... ,7V 
to the notation using g = 5'i,...,5'2 introduced before. In this notation we have 



X, 



(9) 

0-2 



\{g{g + 1) — 2s{s + 1)). In the special case of S' = s we conclude Eq 



1{S{S + 1) - 3s{s + 1)) = -s{s + 1) and, by (50), x[^^ = ~\g{g + 1). Especially 
it follows that x^^^ = 0, which was previously, in section 3, mentioned as a consequence 
of the Wigner-Eckhard theorem. This fact will simplify some calculations in section 5. 
There are two other special cases that should be mentioned: One case is 5 = with 
iV = 1 for integer s. The corresponding eigenstate does not depend on y, see [4], hence 
the spectral curve coincides with its asymptotic line. 

The other case is s = 1/2 and the vector J = (Ji, J2, J3) running through a circle in 
the plane Ji + J2 + J3 = 3j = const with center In this case the eigenvalues of 

H will be constant, i. e. , we have the rare case of isospectrality for spin systems, see [16]. 
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As a further application of the spectral symmetry we calculate the first five moments 
ofH2{y). Let 

TV 

M^{y) = J2{e.{y))\ neN. (55) 

Since the moments can be expressed through the coefficients of the characteristic 
polynomial, see, e. g. , [17] (14.35), the M„(y) will be polynomials in the variable y 
of maximal degree n and their coefficients will exhibit the same reflection symmetry as 
the coefficients of p{x,y), see (29). For example, 

M,{y) = M,o{l + y) . (56) 
The coefficient Miq can be calculated directly or by using 

N 

M,{l-y)= J2{Eo-e.{y)) (57) 

= NE,-M,{y) (58) 

= iVEo-Mio(l + y) (59) 

= Mio(l + l-y), (60) 

hence 

Mio = ^iVEo. (61) 
Similarly, 

M2{y) = M2o(l + y^) + M^iy , (62) 

M2(l -y) = (2M20 + M2i)(l - y) + M20I/' (63) 

N 

= J2{E, - eM? (64) 

= NEl- 2Eo Mi{y) + AUy) (65) 

= NE^o - 2EoMio(l + y) + M2o(l + y^) + M211/ , (66) 

hence 

M2{y) = M2o(l + y') + (^oMio - M2o)y • (67) 

Here and in the following equations we assume the Mno as given since they can be 
expressed through the well-known eigenvalues of H2{0). By the same method we obtain 

Ms{y) = M3o(l + y') + ^{EoM^o - Mso)y{l + y) , (68) 

M^{y) = M4o(l + y") + 2(EoM3o - M^^)y{l + y^) 

+ (-E^Mio + 6E2M20 - 8E0M30 + 3M40) ^/^ (69) 

and 

M^{y) = M5o(l + y') + \{EoM,o - M^,)y\l + y^) 

+ (-^o'^io + 5^0^20 - 5^o'^3o + M50) y"" • (70) 
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For n > Q this method fails to give M„2 in terms of Mfco, k < n and simple formulas 
like those above seem to be no longer available. Note, however, that the first five 
moments have been derived solely by utilizing y = results and spectral symmetry. 
The knowledge about the explicit matrix representation of -^2(1/) has not yet been used. 
The first five moments of H2{y) depend on S only implicitly through Eo and Mno- If 
these arguments are explicitly evaluated it is possible to calculate the corresponding 
total moments of the 2-chain by summation over S and to compare the results with 
those in the literature, see [8] and [9]. This has been done successfully but will not 
be detailed here. Rather, in the next two sections, we will use an alternate method to 
calculate moments which leads to results not already reported in the literature. 



5. Moments of the two-chain 



The total k-th moment of H2{y) will be defined as 
Tr H2(y)'' 

M.) - . (71) 

Obviously, it is a polynomial in the variable y with the same reflection symmetry 
as the local moments Mk{y) defined in (55). The knowledge of the first terms in 
the perturbation series of the eigenvalues eg^siy) of -^2(1/) will give the corresponding 
coefficients of tk{y) via summation over the quantum numbers S and g = Si2- Moreover, 
we will see that there is a close connection between these coefficients and the terms of the 
so-called graph expansion of tk{y)- In the case of the two-chain H2{y) this connection 
is even a 1 : 1-correspondence. 
Let 

eg Ay) = 4'^ + 4'^ y + 4'^ y' + 4'^ y' + --- (72) 

denote the perturbation series of the eigenvalues Cg^siy)- Upon multinomial expansion 
we obtain 

eg,siy)' = 4'^' + kxi'^^'-'^x['^y 

j9)(k-2) (9)2 , (9)(k-l) (g)\ 2 



(g){k-3) {g)3 , , /, T\A9)ik-2) ig) {g) , h.A9)i.k-l) (g) \ 3 



+ ... (73) 
Into this series we insert the known equations (50) and (51): 

xi'^ = l{g{g + l)-2s{s + l)) , (74) 

4'^ = liEo - 4'^) = 1{S{S + 1) - g{g + 1) - s{s + 1)) , (75) 

= ^4'^ . (76) 
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For x'f^ we will use the general result (114) for tridiagonal matrices and the explicit form 
(24) of the secondary diagonal entries of H2{y)- This yields, after some calculations, for 

1 



0^2 



16(2(7+1) 

{g-2s-l){g + 2s + l){g-s-S-l){g + s-S){g-s + S){g + s + S + l] 

^7(2^?-l) 

{g - 2s){g + 2s + 2){g - s - S){g + s-S + l){g -s + S + l){g + S + S + 2) 

{g + l){2g + ?>) 



(77) 

= - -x^^^ + ^ + — (1 + Hs + 1)) + - , ''^\ , (78) 

32 ° 16 256^ ^(^ + 1) ^g{g + 1) - 3 ' ^ ' 

where 

= 1(1 + 2s)\s - S)\\ + s + S)^ (79) 
16 

and 

^2 = ^(1 + 4s)(3 + 4s)(-l + 2s - 2S){\ + 2s - 2S){\ + 2s + 25)(3 + 2s + 2S) .(80) 

For (7 = the first term in the brackets in (77) or, equivalently, the term g(^\x) (^8) 
has to be set to 0. These results have to be inserted into 

2s s+g 

^M = Y. E (25 + l)e.,5(y)^ (81) 

9=0 S=|s-g| 

The double sum in (81) can be evaluated up to third order in y by adopting the following 
strategy: First, expand all terms depending on S[S + 1) into polynomials in and 
express the resulting sums of the form 

(25 + 1)^0-^, (82) 

S=\s-g\ 

up to a factor 2(7 + 1, as polynomials in Xg^^ This makes is possible to write the final 
sums Xlglo • • •^0^'*'*' ill terms of the known moments t\{^^ (dimer polynomials). 

Before writing down the final results of the above-sketched calculations we will 
shortly recapitulate the general theory of moments of Heisenberg Hamiltonians if, 
following [8] and [9], that is also the starting point of the high temperature expansion of 
the corresponding thermodynamic quantities. The /c-th moment of ii has the following 
"graph expansion" 

TriJ^- 

The latter summation extends over a finite set of (multi-)graphs with k edges (bonds). 
Q denotes a certain polynomial of the couphng constants of H that is obtained via 
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summing over all embeddings of Q into the spin system, see [9]. p{G) denotes certain 
other "universal" polynomials of maximal degree k in the variable r = s{s + 1) that are 
independent of the coupling constants of the spin system under consideration. 
In our case of the two-chain H2{y) the graphs needed in the graph expansion of t^iy) 

are either of the form « (two vertices and k edges) or (three vertices and a total 

number oi k = ki + k2 edges) . The cases ^ = « oi Q = can be disregarded since 
then p{Q) = 0, see [8]. Moreover, we have Q = 1 + y'^ in the first case of ^ = « and 

g = + y^^ in the second case of ^ = Hence the coefficients of the y-expansion 

of tj^{y) are exactly the universal polynomials p{Q) of the relevant graphs and will be 
written as p{k) in the first case and p{ki, k — ki), ki < k — ki, in the second one. The 
dimer polynomials p{k) are the well-known moments of the one-chain and will be given 
in explicit form in the Appendix B. The polynomials p{ki,k — ki), ki = 2,3,4,5 can 
be expressed in terms of dimer polynomials p{k') as the result of the above-sketched 
calculations. The results are 
k 

p{2,k-2) = — [{-rY + {k-l)r^p{k-2) + krp{k-l)] , (84) 
and 

p{3, k-3) = — [-3k{2r + k- l)(-r)''-^ + 6k^rp{k - 1) 



144 

/k\ . /k\ 1 

85) 



- 2r{6r + k + l)(^^p{k-2)-6r^(^^p{k-3) 

If ki =4,5 the expressions for p{ki,k — ki) are more complex and will involve sums 
over 0{k) dimer polynomials. It will be economic to write the results in terms of the 
modified dimer functions 

P{i)^p{i)-p^, ^ = 0,1,2,..., (86) 
4r + 1 

although the expressions for p{ki, k — ki) will then no longer be manifestly polynomial. 
The final results are: 

p(4, k -4) = 

kr (-32r{4r - 3)(80r + 21) (| - r)^'~^ - 128(-r)'=(r(-30A; + 80r - 9) - 72)j 



69120(4r 

k (P -3k + 2)r{k {12r^ - 5r + 2) - 36r^ + 39r + 6) 



5fcV ^ (A;- l)A:r(16A:2(3r- 1) - 16A;(3r + 4) - 3(284r + 57)) ^ 

^192" ^ ~ 23040 ^ ~ ^ 



2880 

(A; - 3)(A: - 2)ik - l)kr' jQr' - 3r + 2) 

2880 ^ ^ 

fc-3 / . /o \ -«+fc-3 



P{k - 3) 



+ ^ + 1) (^^^(4^ - 3)(16r + 3)^ - ^) ^ + - 2)(4r + l)^(-r)-^+'=-^ j P(z), 

(87) 
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and 

p(5, k — 5) 
Akr'^{Ar - 3) ( ^ - r 



2,. ^, /3 y ^ (80r + 21)(4A: + 8r-7) 



46080(4r + 1) 



64A;(-r 



20A;r2 - 45fcr - 15fc + IGOr^ - 98r2 - 39r + 15 



46080(4r + 1) 



7,2^/, , A;2(96r-32) + A;(-96r- 68)- 3(532r + lll)^,, 
k^rP{k - 1) + (A; - l)kr '- ^ — '- ^ '-P{k - 2) 



384 ' ' ^ ' 46080 

2A;(8r - 3)r(A;(2r - 1) - 1) - 108r3 + 372r2 - 141r - 45 
-[k- 2)[k - l)k 11520(8r-3) ^^'^ " 

8A;r2 - 5A;r + A; + 24r^ - 52r2 + 34r + 6 
-{k-m- m - ^)kr P{k - 4) 

4r^ — 3r -I- 1 

-{k-A){k- 3){k -2){k- l)kr' ^^^^Q^ P{k - 5) 

k .V 3 „ /3 \-i+k-4 



3840 

i=0 



^(i + 1) ( r(4r - 3)(16r + 3)^(4^; + 8r - 7) - r 

.■_n V 256 8 



(4r + l)2(5A;(r- 1) +r(6r- 17) + 5)(-r)"'+^-3 ]P{i). 



The formulas (84), (85), (87) and (88) have been checked with the known results 
for /c = 2, . . . , 8 in [8] and unpublished results for k = 9,10 that have been obtained by 
different methods. 



6. Moments of the general triangle 

The graph expansion of tk = j^^^^^Ti for the general triangle Hamiltonian (1) 
involves additional triangular graphs with /ci, A;2, k2 bonds such that ki + k2 + k^, = k. 
We will write the corresponding universal polynomials as p{ki, k2, k^), ki < k2 < k^. 
Following the ideas of the last section, a possible strategy to determine these polynomials 
would consist of expanding tk{z) into an even power series w. r. t. z, if the Hamiltonian 
(1) is written in the form 

H{z) = Jg Si ■ S2 + J (s2 ■ S3 + Si ■ S3 j + ( J3 - J) Z (^82 ■ S3 - Si ■ 53^ . (89) 

In this paper we will take only the first step of this expansion and consider the integrable 
case z = 0. Even this case is sufficient to obtain explicit expressions for some classes of 
universal polynomials. 

Consider the eigenvalues Cg^s of -^^(0) that, according to (50) and the remark following 
(25), can be written as 

eg,s = J3xi'^ + 2Jx['^ (90) 
= J,xi'^+J{Eo-x[^^). (91) 
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e'^^=T.(f)j'^^oJ'^'iEo-xl^^r' (92) 



Hence 

k 

= E E 0) C " j'-'Ei;{-ir^-^xi^^^'~^^ . (93) 

1=0 ^,=o \ l^ J 

The moments tfc(O) are obtained by summing (93) over S and (7. This leads to sums of 
the form ^^'^^^_^|(2S'+ l)-E'o • 1^ can be shown that these are of the form (2s + l)(25f + 1) 
times a polynomial in the variables r = s{s + 1) and g{g + 1). Hence we obtain the 
following expansion 

2 {2S+l)E^,={2s + l){2g + l)Y,P,,x'i^\ (94) 

S=\s-g\ A=0 

which implicitly defines the polynomials P^x{r). For our applications we will only use 
the following ones: 

P„„ = 1, n = 0,1,2,..., (95) 

Pio = 0, (96) 

P20 = \r\ P21 = \r, (97) 

Pm = - ^r^, P31 = lr(6r - 1), P32 = 2r . (98) 

Note that the factor (2s + l){2g + 1) is just the degeneracy of the eigenvalue Xq^ of 
Si ■ S2 in the total Hilbert space. Hence the sums of the form Ylf=o (2s + l){2g + 1) x'^^^ 
are equal to (2s + lYp{\) and we obtain 

^^(0) = E E E C 7 ^3^'-'(-l)'-'-^^M P(fc + A - ^) . (99) 

£=0 M=0 A=0 ^ ^ 

We will only consider two terms of the expansion (99), namely those with i = k — 2 and 
i = k — 3. Let [^(x)]^; denote the coefficient of x in the polynomial A{x). Then the first 
term is 

2 M 



[^fc(0)]j-^,2= (-l)'-^P,AP(fc + A-/i) (100) 

^ ^ /t=0 A=0 ^^-^ 

= pil,l,k-2) + 2p{2,k) . (101) 
After some calculations, using (84) and (85), we finally obtain 

p(l, 1, A; - 2) = ^ [-{-rf + {k - l^^k - 2) + {k - 2)p{k - 1)] , 

for A; = 3, 4,... . (102) 
For ^ = k — 3 the analogous results read 

k- 3 , 



[4(0)]4-3,3= (-l)^-^P,AP(A: + A-/i) (103) 

^ ^ ^=0 A=0 ^^'^ 

= 2(p(l,2,A;-3) + 2p(3,A;)) . (104) 
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and 



pil,2,k-3) = ^{{k + 2r- 
+ r [{k-l) ((2- 
- 2kp{k- 1)]} 



l)(-r)^-i 

- k)rp{k - 3) + {-k + 2r + 3) p{k - 2)) 
, for A; = 4, 5, ... . 



(105) 



The formulas (102) and (105) agree with the values of the polynomials published in 
[8] for k = 4, . . . , 8 and with unpublished results for /c = 9, 10 obtained by different 
methods. For an example see Appendix B. 

It should be mentioned that the universal polynomials (84), (85), (87), (88), (102), 
and (105) do not only occur in the graph expansion of the moments of the two-chain or 
the general triangle but for arbitrary Heisenberg spin systems. As a rule, general formu- 
las for these polynomials are rare. The only ones previously known to the author are, 
besides the explicit expressions for the dimer polynomials p{k), see Appendix B, those 
which result from certain rules given in [8]. For example, p{Q) = if ^ is a single-bond 
connection of two disjoint (possibly empty) graphs. Another rule in [8] refers to the 
case where Q contains two vertices connected by a single-bond chain of length n > 2. 
Let be the graph resulting by shortening the single-bond chain by one vertex. Then 
p{G) = ^p{Q), where i is the total number of bonds of Q. (Note that the authors of 
[8] use a definition of the graph polynomials p{Q) that differs from our definition by 
a multinomial factor). In particular, the last rule implies that the polynomial g„ of a 
single-bond polygon with n > 2 vertices will be qn{r) = -^r^r"'. Of course, this rule can 
also be applied to graphs that can ultimately be reduced to those for which (84), (85), 
(87), (88), (102), and (105) apply. 

7. Summary and outlook 

The general quantum spin triangle is not integrable in the strict sense that its eigenvalues 
and eigenvectors can be given in closed form, as, e. g. , for the isosceles triangle. However, 
our results show that we are rather close to integrability. The reduced Hamiltonian has 
a triangular matrix representation the entries of which are explicitly given. Hence also 
the characteristic polynomial can be expressed either recursively or explicitly in terms 
of finite sums, see Appendix A. Similarly, the perturbation series for the eigenvalues 
could, in principle, be calculated up to any finite order, although with rapidly increasing 
complexity. Therefore we have contented ourselves with third order perturbation results. 
Comparing this situation with, e. g. , the Bethe ansatz for the s = | chain with spins 
shows that the difference is not huge. Also in the Bethe ansatz the eigenvalues have to 
be calculated as the numerical solutions of certain equations which become more and 
more complex for large A^. Perhaps the general triangle should also be regarded as 
"quantum integrable" in some weaker sense. 

One virtue of integrability in whatever sense is the possibility to calculate certain limits. 
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e. g. , the thermodynamic hmit — > oo in the case of the s = ^ 



I chain. In the present 



case a possible candidate is the classical limit s — oo, see [18] for the equilateral triangle 
case. This seems to be a sensible project for future work. Another aspect that has only 
be touched in the present paper is the closer investigation of physical properties of the 
general spin triangle, as, for instance, ground state diagrams or magnetization curves. 
As already pointed out in the Introduction, these properties could be of paradigmatic 
character in order to understand more complex systems. 
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Appendix A: Properties of tridiagonal matrices 

We will collect some properties of tridiagonal matrices relevant for this article without 
claiming originality. For sake of simplicity we assume T to be an x A^-symmetric 
real tridiagonal matrix with positive entries in the secondary diagonals. Let T*^"^ denote 
the submatrix of T obtained by deleting its first n rows and colums and D„ = det T^'^\ 
Then Dq = det T satisfies the obvious recursion relation 



It follows that the explicit formula for det T does not contain A^! terms as for general 
T, but only -F/v+i terms, where denotes the A^-th Fibonacci number. For exam- 
ple, if A^ = 3 we have det T = T11T22T33 — — Tf2^33 which consists of F4 = 3 
terms instead of 3! = 6 ones. The number of terms still grows exponentially with 
A^ but not super-exponential like A^!. Consequently, the computer-algebraic calculation 
of the characteristic polynomial of T is markedly faster than for a general A^ x A^-matrix. 

We will write down the explicit formula for det T in order to substantiate our claim 
that the characteristic polynomial of the general triangle can be given in closed form. 
To this end let {di, . . . , cIn) denote the diagonal of T and (c2, . . . , cat) its secondary 
diagonals. Further let B{m, k) denote the set of (0, l)-sequences of length m + k 
containing exactly m zeroes and k ones. There are {"^^^) such sequences. For any 
a G B{m, k) we define 



(106) 




(107) 
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Then 
det T 



m+k 



TO, fc e N 
m + 2k^ N 



E n 

aGB{m,k) i=l 



d, 



.2 

''n{i,a) 



ait) 



0, 
1. 



(108) 



Another simphfication occurs for the perturbation series of, say, a diagonal matrix 
A with pairwise different eigenvalues perturbed by a symmetric tridiagonal matrix y B. 
We will write the tridiagonal matrix A + y B in the form 



( * 


* 








* 


* 


* 








* 


* 


yy/Bn-2 








yVBn-2 


An-i + ya„ 











yy/Bn-l 






















































yVBn-1 

An + ya^ 

yy/Bn+l 







\ 







yVK:;'! 

An+1 + ya^+i * 

* * * 

* * * 

* * y 



(109) 



which is symmetric w. r. t. the n-th diagonal entry (note that there is no Bn)- This 
symmetric form facilitates the notation in what follows. We introduce, for fixed n, the 
abbreviations 



ai = Ai 

A - 



k -= 

Let 



Aj — A^ 
Bi 



Ai Ar. 



[i ^ n) 
[i ^ n) 



any 



^y" 



(110) 

(111) 

(112) 



(113) 



denote the perturbation series for the n-th eigenvalue x^"^ oi A + yB, where we have 
anticipated the elementary result for the first two terms. Then the next terms are of 
the form 



X 



X 



in) 
2 

(n) 



X 



in) 



- K-i 

bn-2bn-l 



Jn+l 



-1 + RT, 



+ 



ft2 



bn_if3l^, + RT . 



(114) 
(115) 

(116) 



fln-l fln-l fln-l 

Here "RT" (reflected terms) stands for the same terms as before but with reflected in- 
dices n — i -v^ n + i. The above formulas also hold for small n = 1, 2, 3, . . . or large 
n = N,N — 1,N — 2, . . . if the corresponding A^ and Bi are set to zero. It is remarkable 
that the perturbation has a "local" character in so far as x^^J and x^J_^_-^ only depend 
on matrix elements with indices n — k, . . . ,n + k. Therefore it is possible to obtain the 
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above results by computer-algebraic calculations with at most 5x5- matrices. We note 
that there are further simplifications for our case of A + yB = si ■ S2 + y S2 ■ S3, see 
(50)-(54), which do not hold in the case of general tridiagonal matrices. 

A further property of tridiagonal matrices T with non-vanishing entries in the 
secondary diagonals is that their eigenvalues are never degenerate. 
For the proof, let (p be an eigenvector of T with eigenvalue x, that is 

Aip = {T-xl)ip = . (117) 

It follows that Aiicpi + Ai2ip2 = 0, or ip2 = since A12 = T12 7^ 0. Similarly, 

+ ^22V^2 + ^23V'3 = implies 

<^3 = . \ (^11^22 - A12A21) (^1 , (118) 

^12^23 

since ^23 7^ 0. By induction, one easily shows that 

¥^n=(^nAmj (-l)"-^/^("-^Vi, n = 2,...,7V, (119) 

where D*-""^ denotes the n-th principal minor of A, e. g. , D^^^ = An, D^"^^ = 
A11A22 — A12A21, and so on. (119) implies that the eigenspace of T corresponding 
to the eigenvalue x will be one-dimensional, which completes the proof. 
As a by-product we have obtained an explicit formula for the eigenvector ip corresponding 
to the eigenvalue x of a tridiagonal matrix T. 



Appendix B: Moments of the one-chain 



Since the eigenvalues of the one-chain (or "spin dimer") are 

cs = \{S{S + 1) — 2s(s -|- 1)) with degeneracy 25* -|- 1, its k-th moment is given by 



p{k) 



(2s + 1) 



1 

^ 5=0 



2s(s + l)) 



(120) 



This sum can be conveniently evaluated by computer- algebra software but it is not an 
explicit formula for the polynomials p{k){r), r = s{s + 1) in the strict sense. However, 
such an explicit formula exists since (120) can be reduced to power sums, which in turn 
are expressible by means of the Bernoulli numbers Bn- 

p{k){r) = 



-3k-2 



k + 



For example. 



k + 1 

V 



k-u+l 2{k-u+l) 



J 

2/i 



(121) 

2{k - Z/ + 1) 
3 



p{2) 



(122) 



p{l, 2, 3) = - {-r\5 + 2r) + r [5(-4rp(3) + (-3 + 2r)p(4)) - 12p(5)]} 
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p(3) = - (123) 
o 

p(4) = -3-r^(2-2r + 3r2) (124) 

1 

p(5) = - -r\2 -2r + 2r^) . (125) 
6 

From these polynomials we can, for example, calculate p(l,2,3) by using (105): 

(126) 

= V(l - 2r) , (127) 
in accordance with [8] and [9]. 

Let us write p{k){r) = X]n=2 Then the first few leading coefficients of the 
polynomials p{k) assume the following form: 

f ^ : , even, 

] : /c odd, ^ ' 

{k(2-k) , 
: A; odd, ^ ' 

( (k-2fK7k+2) . , 

J.k) _ I 720{fc+i) • «;even, 

"fc-2 — \ (k-3)(k-l)(7k~5) , T n K^'^^J 



(.-3)(fe-l)(7.-5) . j^^^^ 
(4-fc)(/t-2)fc(62fc3_l74fc2+fc+l74) 



k even, 



(fc) ^ J 60480(fc+l) • «ven, 

*^-3 j {3-A:)(fc-l)(62fc3-360fc2+523fc-105) ^ 



60480 



odd. 



Another representation of the polynomials p{k) can be obtained from the Euler- 
MacLaurin summation formula, see [19] 23.1.30. After some simplifications we arrive at 
the following expression involving only a single summation: 

^^^^ = {2sliykli 

k 



1 r,l o_2 



8s^ 



.m=l \ / \ \ 

-s(s + l))'=-"'+^F(-m, -A; + m-l ^ 



2' 8s(s + l) 

+ ^ (s^'^ (4A;s + k + {2s + If) + (-s(s + 1))^= (A; + {2s + 1)^)) 



(132) 



Note that the hypergeometric function -F(a, 6; c; z) reduces to a polynomial in 2; if a 
or h are negative integers, see [19] 15.4.1. Hence (132) is a priori a rational function 
of s. Only after the summation it turns out that the factor (2s + 1)^ cancels and the 
remaining sum is a polynomial in r = s{s + 1). 

The p{k) can be viewed as polynomials in the complex variable r. Interestingly, the 
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Figure 5. The sequence of complex zeroes of the dimer polynomials p{k), k = 2,..., 50, where 
different colors correspond to different values of k. Obviously, besides the trivial value of r = 0, there 
are accumulation points for the physical values r ~ s{s + 1), s = 1/2, 1, 3/2, 2, 5/2, 3. 



physical values of r = s{s + 1), s = |, 1, |,2, . . . can be recovered as accumulation 
points of the sequence of complex zeroes of all p{k), k = 2,3,4, .. .. This is not rigorously 
proven but demonstrated by numerical evidence, see figure 5. 
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